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There are various routes for deriving partial radial distribution functions of disordered systems from experi- 
mental diffraction (and/or EXAFS) data. Due to limitations and errors of experimental data, as well as to imper- 
fections of the evaluation procedures, it is of primary importance to confirm that the end result (partial radial 
distribution functions) and the primary information (diffraction data) are consistent with each other. We intro- 
duce a simple approach, based on Reverse Monte Carlo modelling, that is capable of assessing this dilemma. As 
a demonstration, we use the most frequently cited set of "experimental" partial radial distribution functions on 
liquid water and investigate whether the 3 partials (O-O, O-H and H-H) are consistent with the total structure 
factor of pure liquid D2O from neutron diffraction and that of H2O from X-ray diffraction. We find that while 
neutron diffraction on heavy water is in full agreement with all the 3 partials, the addition of X-ray diffraction 
data clearly shows problems with the 0-0 partial radial distribution function. We suggest that the approach 
introduced here may also be used to establish whether partial radial distribution functions obtained from sta- 
tistical theories of the liquid state are consistent with the measured structure factors. 
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1. Introduction 

Experimentally determining the structure of multicomponent disordered materials, such as liquid 
water, aqueous solutions, etc., poses nearly prohibitive difficulties already at the two-particle level, due 
to the lack of measurable information. For instance, in two-component liquids, like water (composing 
atoms: and H) or carbon tetrachloride (composing atoms: C and CI), there are 3 partial distribution 
functions (PRDF), g,;(r): goofr), gOH(r) and g H H(r) for H 2 and gcc(r), gcci(r) and gcicito for CC1 4 . PRDF-s 
cannot be measured separately (one-by-one) in general; instead, composite functions, called total scatter- 
ing structure factors |1] (TSSF), F(Q), may be determined in the reciprocal space. The real-space equiv- 
alent of F(Q) is its Fourier-transform pair, called the total radial distribution function (TRDF), G(r). F(Q) 
and G(r) contain contributions from PRDF-s which depend on the scattering power and concentration of 
the components 

G(r)= £ acjbibjlgijir)-!], (1.1) 
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F(Q) = Pof 47ir 2 G(r) (^p) dr. (1.2) 

o 

In the above equations, c,- and b{ are the molar ratio and the neutron scattering length of species 

1, G(r) is the total radial distribution function, po is the number density and Q is the scattering variable 
(proportional to the scattering angle); indexes i and j run through nuclear species. We would like to stress 
that PRDF-s cannot be considered as primary experimental information; it has been known for some 
time that more than one set of PRDF-s may be equally consistent with a given set of measured structure 
factors Q], containing different isotopic ratios of 1 H and D ( 2 H) in water. Therefore, any particular set of 
the so-called "experimental" PRDF-s should be taken only as one possible interpretation of the measured 
data. That is why the underlying problem of consistency between TSSF-s and PRDF-s deserves particular 
attention. 

Strictly speaking for an unequivocal determination of all the three PRDF-s of any two-component sys- 
tem, three independent measurements would be necessary on the same (at least chemically, if not isotopi- 
cally, identical) material. This requirement makes the task difficult and performable only for a relatively 
small number of systems. For this reason, it is of utmost importance to be aware of the quality/reliability 
of the PRDF-s derived. 

Neutron diffraction with isotopic substitution |4, 5] (NDIS) and anomalous X-ray scattering al- 
low, in fortunate cases, to separate all the PRDF-s. These fortunate cases are all two-component systems. 
Arguably the best known example of isotopic substitution is replacing H by 2 H (deuterium). The trick 
was tried first on water |6], which liquid has been since probed many times by NDIS (see, e.g., 
The fact that one single experiment has not been able to clarify water PRDF-s for once and forever in- 
dicates that there must be some unresolved issues concerning the procedure. Indeed, a most (in)famous 
property of 1 H nuclei (protons) is their incredibly large incoherent inelastic scattering cross section for 
neutrons 1 9]. This property renders more than 95% of the signal measured by neutron diffraction useless 
("background-to-be-subtracted") from the structural point of view (for a more detailed exposition of the 
problem, see (3J]). 

One of the latest compilation or PRDF-s resulting from H/D NDIS experiments on water can be found 
in Jioh : it is rather helpful that numerical data for PRDF-s resulting from this work are posted also on the 
Internet Hill . Due to difficulties mentioned above it would be rather comforting if an independent (and 
preferably, a positive) assessment of these results could be provided. 

In this work, we propose a possible approach. As primary data, recent neutron diffraction results 
from the same author |8] on liquid heavy water, as well as the most recent published set of X-ray data of 
Fu et al. (l^] will be applied. The procedure uses the so-called Reverse Monte Carlo (RMC) algorithm Hill , 
based on which one of the present authors (LP) has published, together with Prof. Orest Pizio, a related 
approach flUl . that later could be applied for assessing the performance of interaction potentials of liquid 
water Hill . 

2. Reverse Monte Carlo modelling 

Reverse Monte Carlo jl5II3Jl3l is a simple tool for constructing large, three-dimensional structural 
models that are consistent with total scattering structure factors (within the estimated level of their er- 
rors) obtained from diffraction experiments. Via random movements of particles, the difference between 
experimental and model total structure factors (calculated similarly to the j 2 -statistics) is minimized. As 
a result, by the end of the calculation a particle configuration is available that is consistent with the ex- 
perimental structure factor(s). If the structure is analyzed further, partial radial distribution functions, 
as well as other structural characteristics (neighbour distributions, cosine distribution of bond angles) 
can be calculated from the particle configurations. 

A possible algorithm that can realize the above features may be outlined as follows H3l : 

1. Start with an initial collection of particle coordinates in a cubic box; this may be a crystalline or a 
random distribution of-at least a few thousands of-particles, or even the final particle distribution 
from a previous simulation. 
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2. Calculate the partial radial distribution functions for the configuration. Compose total radial distri- 
bution functions, according to the experimental weighting factors. Use Fourier transformation for 
calculating total scattering structure factors. 

3. Calculate differences between model and experimental functions as follows (shown here for one 
single TSSF): 

i ° 2 

The "C" and "E" superscripts refer to "calculated" and "experimental" functions, respectively; o is 
a control parameter that is related to the assumed level of experimental errors. 

4. Move one particle at random. 

5. Calculate PRDF-s, TRDF-s, TSSF-s and from them, also the x 2 , for the new position. 

6. If the x 2 for the new position is smaller than it was for the old position (i.e., the difference between 
simulated and measured TSSF-s has become smaller) then accept the move immediately. Otherwise 
accept the move only with a probability that is proportional to exp(-Aj 2 ); accepting "bad" moves 
with such small but finite probability will prevent calculations from sticking in local minima. If a 
move is "accepted" then the "new" position becomes the "old" one for the next attempted move. 

7. Continue from step 4. 

The most valuable feature of the RMC method concerning the present purposes is that it can incor- 
porate any piece of information that can be calculated directly from particle coordinates. Partial radial 
distribution functions from traditional evaluation of diffraction data | 4, 5] fall into this category. In this 
case, if consistency with all input data is reached then it may be stated that these input data are mutually 
consistent. If, however, some of the input data cannot be approached within their uncertainties then it 
means that parts of the input data set are not consistent with other pieces of input information. In our 
case this would mean that some of the input PRDF-s from the literature 1 10] would not be consistent with 
the experimental input total scattering structure factor(s). 

In the RMC calculations that are the basis of the present work, one total scattering structure factor 
from neutron diffraction 121131, one TSSF from X-ray diffraction [2M an d three partial radial distribution 
functions from [HI are ap plie d as input data for each calculation. The PRDF-s appear for various 
thermodynamic conditions in liollllll : here the ones denoted as set "(a)" are being used. Data for set 
"(a)" were obtained at 298 K and 0.1 MPa (ambient conditions); the molecular number density was 0.0334 
molecules A -3 . 

We report two calculations: as real experimental results, the first one uses neutron diffraction data 
only whereas in the other one, neutron and X-ray TSSF-s are both applied. Naturally, both calculations 
involve the same set of PRDF-s. In both RMC computations, the simulation box contained 2000 water 
molecules (6000 atoms). Goodness-of-fit values, R w -s (which are basically sums of the squared differences, 
see below), are reported in a normalized form, so that variations in terms of the number of r and Q 
points considered would not affect the assessment; additionally, the applied r and Q ranges were kept as 
uniform as possible. 

For making the definitions of the different R w -s used throughout this study clearer, below we provide 
the appropriate equations as to how to calculate these differences (for PRDF-s, only the example of the 
0-0 g(r) is shown): 

2 _ £,(f c (Q,)-f e (Q,)) 2 

i?w [F(Q)] - L . F E (Q . )2 • W 

R W [goo(D] - . (2.3) 

In the above expressions, AT,- and Nj are the number of Q and r points, respectively, for the experi- 
mental TSSF and "experimental" g(r)-s, respectively. Indices "C" and "E" refer to "RMC calculated" and 
"experimental" quantities. 



23606-3 



Z. Steinczinger, L. Pusztai 



3. Results and discussion 

Table [l] summarizes goodness-of-fit values for each calculation reported here. Just by looking at the 
numbers, it becomes evident that the scheme proposed here looks rather promising by being decisive in 
a couple of important issues. It is clear that consistently it is the O-H PRDF that is the least consistent 
with primary experimental data (the TSSF-s); this is in line with previous findings am. 



Table 1. Goodness-of-fit (R w ) values for individual data sets [F^Q), F^-(Q) and the three partial g(r)'s]. 
"R w sum" is simply the sum of individual R w values. (ND: neutron diffraction; XRD: X-ray diffraction.) 





ND only 


ND+XRD 


R w [ND F(Q)] 


4.67% 


4.67% 


R w [XRD F(Q)] 




4.27% 


Rw [goo(r)] 


4.03% 


8.79% 


Rw [gOH(r)] 


16.91% 


16.93% 


Rw [gHH(r)l 


3.87% 


3.90% 


R w sum 




38.56% 


R w sum (without XRD) 


29.41% 


34.29% 



It is instructive to look at figures[T]and[2] so that the level of consistency between experimental TSSF-s 
and PRDF-s may be visualized better. We note here that no separate figures are provided for calculations 
where only neutron diffraction data were considered, since the tiny differences in terms of the goodness- 
of-fits are not visible, apart from the 0-0 PRDF, which is discussed below more in detail. 




Figure 1. RMC modelling "experimental" partial 
radial distribution functions of liquid water (set 
"(a)" in (l^, [l3l)> together with neutron diffrac- 
tion TSSF of heavy water 0H] (part a) and X-ray 
diffraction TSSF of light water (part b). Solid 
line: experiment; dotted line: RMC. 



Figure 2. RMC modelling "experimental" O-H (part 
a) and H-H (part b) partial radial distribution func- 
tions of liquid water (set (a) in (l(l[ll]]), together 
with neutron diffraction TSSF of heavy water (tHU 
and X-ray diffraction TSSF of light water (ill. Solid 
line: "experiment" ( (ll|l); dotted line: RMC. 



Even though in each calculation the R w factor for the O-H PRDF is significantly higher than that for 
the other two partials, agreement between RMC and "experiment", and consequently, between primary 
experimental information and PRDF-s derived from it in [10], is rather convincing. 

It is obvious from both tableland figure[3]that the introduction of X-ray diffraction data influences 
the 0-0 PRDF the most. Clearly, the quantitative consistency with the neutron data (i.e., where the calcu- 
lation did NOT contain X-ray data), has worsened substantially. This finding indicates that PRDF-s derived 
from neutron diffraction data with isotopic substitution are not entirely consistent with the latest (and 
independent) set of X-ray TSSF. Note also that the X-ray data are dominated by the 0-0 PRDF (whereas to 
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the neutron weighted TSSF of D 2 0, the contribution is less than 10 %, see HQ, SHU), so triat -P rov i ded 
that errors in the two TSSF-s are comparable-in terms of 0-0 correlations, the X-ray TSSF should be 
considered as decisive. This is a straightforward indication that further improvements would be needed 
to the procedures applied during data evaluation of neutron diffraction data sets that work with H/D 
isotopic substitution. 

Given the simplicity of the approach, its application might be fruitful in other areas, too. For example, 
statistical theories have also been frequently used for producing PRDF-s of various liquids (see, e.g. (ljj 
Eoll and references therein). To assess the consistency level of these theories with experimental data, the 
RMC-based calculations described above might be adequate. 



4. Conclusions 

A simple approach has been devised in order to establish the level of consistency between experimen- 
tal diffraction data, that can be obtained in the reciprocal space in the form of composite total scattering 
structure factors, and partial radial distribution functions that may be derived from them. 

For the particular case of the most frequently cited liquid water PRDF-s (lo|] it has been established 
that they are consistent with the most reliable neutron diffraction data set, taken on pure heavy water 
tilla]. On the other hand, the 0-0 partial radial distribution function provided by 1 10] is not entirely in 
agreement with the results of the latest X-ray diffraction measurements (l2|] on liquid water. 
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He3a/ie>KHi/m, 3arasibHi/ift MeTOfl fl/ia nepeBipKi/i y3rofl>i<eHHfl 
Ai/i<|>pai<LJiiMHiix flaHMX3 napi^ia/ibHMMM paflia/ibHi/iMi/i 
<t>yHi<LVflMi/i po3nofli/iy, ai<i OTpi/iMam 3 hi/ix: npi/ii</iafl pifli<oT 

BOflM 

X. LUTawH^HrepP, J]. llycTaP 

LLlKO/ia flpyroro oyneHfl Eyflai Haflb Aman, H-1 121, EyflaneuiT, yropmnHa 

iHCTHTyr c|)i3HKW i onTWKH TBepfloro Ti/ia,BifHepiBCbKnti floaiiflHufi uemp (t>i3WKH, YropcbKa AKafleMm HayK, 
H-1 525 ByflaneujT, YropmwHa 

ICHye fleKi^bKa ui/wxiB OTpuMaHHfl napuja/ibHux paflia/ibHwx cpyHKi4iPi po3nofli/iy b HeBnopaflKOBaHi/ix cucre- 
Max i3 eKcnepwMeHTa/ibHwx flwct>paKUMHi<ix (i/a6o EXAFS) flaHwx. Hepe3 o6Me>KeHiCTb Ta noxn6KH eKcnepi/iMeH- 
Ta/ibHux flaHnx, hk i HeflOCKOHa^icTb o64ncmoBa.nbHnx npoueflyp, nepuioneproBOi' Bax/inBOCTi Ha6yBae nifl- 
TBepflxeHHfl toto, luo KiHi4&Bi pe3y/ibTa™ (napuja/ibHi paflia/ibHi cbyHKi^iT po3nofli/iy) i nepBUHHa mefiopMaujfl 
(flncppaKi4ifiHi flam) y3roflxyK)TbCfl Mix co6ok>. nponoHyeTbca npocTi/m niflxifl, hk\/\v\ 6a3yeTb«i Ha o6epHeHOMy 
MOfleyiiOBaHHi MoHTe Kap/io, akhw cnpoMOXHufi po3B'fl3a™ u,k> flu/ieMy. B HKOcri fleMOHCrpauji mh BUKopucTO- 
BycMO Ha tiLinTOBa Hi lu n^l Ha6ip "eKcnepuMeHTa^bHux" napuja/ibHux paflia/ibHux cpyHKL^iiTi po3nofli/iy fl/w BOflu 
i floaiiflxyeMO, hh yci Tpn napuja/ibHi po3nofli/iw (O-O, O-H i H-H) y3roflxyK)TbCfl 3 noBHUM CTpyKTypHUM cf>a- 
KTOpoM 4mctoT BOflu H2O, OTpi/iMaHWM i3 fli/icbpaKi^iY X-npoMeHiB. Mn noKa3yeMO, mo X04a flaHi HePiTpoHHOro 
po3ciK)BaHHfl Ha BaxKifi BOfli uJ/ikom BiflnoBiflaKDTb yciM napuja/ibHUM po3nofli/iaM, flOflaTKOBe BpaxyBaHHfl 
flaHHX po3ciflHHfl X-npoMeHiB BWRBnne npo6/ieMW 3 napuja/ibHOK) cfiyHKujeK) po3nofli/iy O-O. Mn nponoHye- 
mo 3acTOCOByBaTH 3anponoHOBaHnPi TyT niflxifl TaKOx fl/ia BWHCHeHHH toto, 4H napuja/ibHi paflia/ibHi cpyHKi^iY 

p03nOfli/iy, flKi OTpHMaHi i3 CTaTHCTWHHUX Teopifi piflKOI"0 CTaHy, y3l"OflXyK)TbCfl i3 BWMipflHUMW CTpyKTypHUM n 

cfiaKTopaMH. 

ICmowoBi cnoBa: HempoHHe P03cirhhr, napujanbHi paflia/ibHi cpyHKLfif po3nofiijiy, o6epHene MOflejuoBaHHH 
MoHTe Kapjio 
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